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ABSTRACT 

We analyze non-Gaussianity (NG) due to the primordial bispectrum and trispectrum using 
cosmic microwave background temperature maps of WMAP 7-year data. We first apply the 
perturb ativc formulae of Minkowski functionals (MFs) up to second-order NG derived by 
iMatsubar a (2010i). which enable us to give limits on cubic NG parametrized with tnl and 
^NL as well as various types of quadratic NG parametrized with /nl- We find no signature 
of primordial NG in WMAP 7-year data, but give constraints on the local-type, equilateral- 
type, orthogonal-type /nl: f^l''> = 20 ± 42, /^^^^ = -121 ± 208, and fl,f'> = -129 ± 171, 
respectively, and tnl/IO** = —7.6 ±8.7, and (7nl/10^ = —1.9 ±6.4. We also find that these 
constraints are consistent with the limits from skewness and kurtosis parameters which 
characterize the perturbative corrections of MFs. 

Key words: Cosmology: early Universe - cosmic microwave background - methods: 
statistical - analytical 
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i INTRODUCTION 

■Non-Gaussianity (NG) in primordial curvature perturba- 
tion is a key observational probe to study physics in 
parly Universe. Single slowly-rolling scalar inflation model 
predi cts 



too 



199C; 



120031 : 



Falk et al 



small 
r; 



N G to be observed ( Salopek fc Bom 



1993"; 'Ga ngui et all 1 1994 

Maldacena i2003i : iBartolo. Matarrese 
■Variants of the simple inflationary mo dels generate much 
}iighe r levels of NG: multiple fields 



Acquaviva et al.l 
fc Riottd |200^ 



. (|Linde fc Mukhanovl 

Il997l : iLvth Ungarelli fc W ands 200 3); niodulat ed reheating 
scenario s (iDvali. Gruzinoy fc Zaldarriagal 120041') : 



flation (iGupta et al.l |2002|: IMoss fc Xior^ 



warm m- 



. 1200711: ek p vrotic 

model (iKovama et all l2007l: ICreminelli. fc Senatord l2007l : 



iBuchbinder. Khourv fc OvrutI bOOTl V These NGs have been 
constrained using the cubic-order statist ics (i.e., bispectrum 



from CMB temperature maps of WMAP (iKomatsu et al.|[201 



lOkamoto fclb3l2002l ). So-called local- type NG is described as 
the nonlinear correction in Bardeen's curvature as follows: 

<E> = + /nl(0'-((^'))+W'-, (1) 

(e.g., iKogo fc Komatsul |2006| ) where (p is an auxiliary Gaus- 
sian variable and $ is related to the curvature perturba- 
tion ^ as $ = 3/5C. In this local-type m odel, tnl is 
equal to 36/nl/25. ISuvama fc Yamaguchil (|2008l ) find that all 
classes of multi-inflation models satisfy the inequality condi- 
tion Tnl > 36/nl/25 derived from the Cauchy-Schwarz in- 
equality. The relation between tnl and /nl is therefore an- 
other powerful probe of inflation models (e.g., ISuvama et al.l 
I2OI0I ). Current WMAP constraints on tnl and (/nl are ob- 
tained by several groups: -5.6 < g^h/W^ < 6.4 (95%C.L.) 
from the N-point probability density function jVielva fc Sand 
-12.4 < anl/IO^ < 15.6 (9 5%C.L.) using trispectrum 



Creminelli. Senatore. fc Zaldarriagal |2007| : 
20081 : ICurto et al.ll201ll )" 



Yadav fc Wandeltl 



3 Ml: 



So far NG analysis has been mainly focused on the mea- 
surement of primordial bispectrum. Primordial trispectrum, 
the next higher-order term, also provides an important probe 
to differentiate inflation models. Trispectrum is generally de- 
scribed into two different kinds of connected parts and their 
amplitudes is commonly parametrized by (;nl and tnl (e.g. 



hikageOkmi . nagoya- u . ac . j p 



( Fergusson. Regan. Shellaxd! l2010l h -0.6 < tnl/IO^ < 3.3 
(95%C.L.) and -7.4 < anl/IO^ < 8.2 ( 95%C.L.) using kur- 
tosis power spectrum (ISmidt et al.l 2010lh Pla nck is expected 
to reach Atnl ~ 560 (jKogo fc Komatsull2006l '). 

Full analysis of trispectrum is computationally very ex- 
pensive because their number of configurations are enormous. 
Minkowski functionals (MFs) are another powerful tool to con- 
strain NGs from the aspect of morphological properties of the 
density structure. MFs and polyspectrum such as bispectrum 
and trispectrum are very different statistics. MFs are real-space 
statistics, whereas the polyspectrum is defined in harmonic 
space. MFs naturally incorporate the information of all orders 
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of polyspectra and thus they are complementary to the com- 
monly used measurement using polyspectr a. Accordin g to th e 
perturbative formulae of MFs derived bv iMatsubaral (|2010l ). 
the higher-order NG originated from primordial trispectrum 
is characterized with four kurtosis parameters, which is the 
summation over trispectrum with different weights on tetra- 
hedron configuration. The computation of kurtosis parameters 
is much less expensive compared to the trispectrum estimator. 
MFs are model-independent statistics and hence they have a 
potential to serendipitously find NGs due to unknown sources 
that has not been explored. Consistency check using different 
estimators from the standard one is thereby important to ob- 
tain more robust results. 

In this paper, we first apply the perturbative formulae of 
MFs upto 2nd-order derived by[Matsubara (20lJ) to WMAP 
7-year data of CMB temperature anisotropics. We give new 
constraints on tnl and (jnl as well as local, equilateral, and 
orthogonal types of /nls. This paper is the extension of the 
previous work to give limits on th e local-type /nl from MFs 
of W MAP 3-year datalHikage et al . ( 2008) and BOOMERanG 
data lNatoh et al.l (|201(]f K 

The paper is organized as follows: i n section [21 we re - 
view the perturbative formulae derived bv IMatsubaral (|2010l ). 
In section O we compare the perturbative predictions with nu- 
merical simulations. We apply the perturbative formulae to the 
WMAP data and constrain different NG parameters in section 
m Section [S] is devoted to the summary and discussions. 



2 PERTURBATIVE FORMULAE OF 

MINKOWSKI FUNCTIONALS FOR CMB 
TEMPERATURE MAPS WITH PRIMORDIAL 
NG 

2.1 Minkowski Functionals for Gaussian Fields 

Minkowski functionals (MFs) have been u sed to characterize 
the morphology of a given density field (|Mecke et al.l 1 19941 : 
ISchmalzing fc BuchertI Il997l ) . In 2-dimensional field such as 
CMB temperature maps, three MFs are defined: area frac- 
tion (Vo), circumference (Vi), and Euler characteristics (V2). 
We measure them for the area whose temperature anisotropy 
/ = AT/T normalized by the standard deviation ao = {f'^}^^'^ 
is larger than a threshold value ly. In Gaussian temperature 

(G) 



maps, the fc-th MF is given by 



(2) 



where Hk (u) represent the fc-th Hermite polynomials. The am- 
plitude Ak is given by 

fe 



1 



bJ2 



(3) 



where i^k = 7r'""''^/r(fc/2 -I- 1). The standard deviation ao and 
that of the first derivative ai = (|V/|'^)^''^ are written as a 
sum of the power spectrum C; : 



^^(2/4-1) [l{l + l)YCiW?, 



(4) 



where Wi represents the smoothing kernel determined by the 
pixel and beam window functions and any additional smooth- 
ing. Here we use Gaussian kernel Wi — exp(— + 1)^^/2) and 



9 denotes the smoothing angular scale. Table [T] lists HEALPix 
pixel number A'^ido (the total pixel number is 12N^i^^) and the 
maximum multipole number Zmax at different values of 9. We 
choose the values of A'aide and Zmax so that the effects of pixel 
window and high frequency cut be small. As discussed later, 
measuring MFs of CMB maps with different smoothing scales 
is important to extract configuration dependence contained in 
primordial bispectrum and trispectrum. 



2.2 Perturbative Formulae in Weakly Non-Gaussian 
Fields 



IMatsubaral (120031 ) has applied multivariate Edgeworth expan- 
sion theo rem to deriv e pertu rbative formulae of MFs for a gen- 
eral field. IMatsubaral (|2010t ) extended his analysis and derived 
2nd-order correction of MFs on CMB temperature maps. Ac- 
cording to their works, the MFs of a weakly NG field (i.e. 
ctq <C 1) are written upto the 2nd-order term of ao as 



/2, 



Avkii^) = vl^\iy)ao + vlf' {u)ao 



(2)/ 



(5) 
(6) 



In the following subsection, we review the perturbative formu- 
lae for a CMB temperature map. 



2.2.1 First- order perturbation 

The Ist-order perturbation terms of MFs are characterized by 
three skewness parameters: 

= ^Hk+2H - ^Hu{v) - ^Hk-2{v), (7) 

with 



if) 



Si: 



, Su 



2{|V/pVV) 



(8) 



"0 '^O"! "1 

The skewness parameters are written as the sum of the bispec- 
trum with different weights of triangle configurations: 



Sa = 4 ^ Ihl2l3SAhl2l3Bl,l2l3Wl^Wl2Wl3 



Aiia\ 
where 

Su i^i^i^ 



{h} + {h} + {h} 

^[{hf+{l2}'+{l3f 



(9) 

(10) 
(11) 



12g 



^2{{h}{h} + {hUh} + {h}{h})], (12) 
with q = a\l\f2ao and {?} = l{l + 1), and 



(2/1 + \)(2h + 1)(2;3 + 1) / h h h 







(13) 



Note that the proportional factor of the skewne ss parame- 
ters a re different from the parameters g^'^ use d in IMatsubaral 
l|2003l '): [Hikage. Komatsu fc Matsubaral (|2006l ) and their rela- 
tions are S = S^°\Si = -45'('V3 and Su = -2S^^^ /3. The 
three skewness parameters have different weights of the bispec- 
trum and hence the MFs can extract more information on the 
configuration dependence of bispectrum than using only one 
skewness value of S. The bispectrum Bi-^i^i^ is defined as 
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a;, 



I «i2 



h 

mi 



h 

1712 



h 

ma 



(14) 



where aim is the harmonic coefficients of a given temperature 
anisotropy map. The relation to the reduce bispectrum bi^^i^i^ 



is B,- 



We consider three different types of 



NGs due to primordial bispectra: local type, equilateral type, 
and orthogonal type. Single and multi-field infiation models 
predict local- type NG (eq.[T]), which generates the follo wing 
form of CMB bispectrum (e.g.. iKomatsu fc Spergell200ll 'l: 



B11I2I3 
with 



2f^r^Ii,i2h / r'dr[ai,{r)Pi,ir)Pi,ir) + cyclic], (15) 



(16) 
(17) 



- / k'^dkgTi{k)ji{kr) 
n 



k^ dkP^{k)gTi{k)ji{kr) 



where qti is the radiation transfer function and ji is the spher- 
ical Bessel function. We rewrite /nl as f'^'^^ ■ The local-type 
NG is sensitive to the bispectrum with squeezed configuration 
of triangle wavevectors (Zi <^l2 — l-i). 

Other inflat i on sce n arios with non-canonical kine tic terms 
(|Seerv fc LidsevI l2005l: IChen. Easther fc Liinl bOOTp . Dirac- 
Born- Infeld models (^Alishahiha c t al. 2004P . and Ghost infla- 
tion (jArkani-Hamed ct al. ^04) predict large NG signals in 
equilateral conflguration triangles [li ~ I2 — t-i,)- Equilateral- 
type NG is characterized w ith f'^^^ defined as the amplitude 
of th e following bispectrum (|Babich. Creminelli fc Zaldarriagal 
I2OO4I ): 



where b'^'^-' is a constant value. 

Fig-IDplots the smoothing scale dependence of three skew- 
ness parameters (eq.[8]) from the local-type, equilateral- type, 
orthogonal-type primordial NG components with the unity 
value of /nl and from unmasked point source NG with the 
constant bispectrum b'^'^-' = 10~^^. We add noise and beam 
functions of WMAP 7-year V-|-W co-added maps which be- 
come important at smaller smoothing scale such as S < 10 
arcmin. Scale dependences of skewness parameters are quite 
different among different NG types. Measuring MFs of CMB 
maps with different smoothing scales is important to break 
degeneracy of different NG sources. 



2.2.2 Second-order perturbation 

iMatsubaral (|2010l ) has derived the second-order corrections of 
MFs, which are characterized with the product of skewness 
parameters and four kurtosis parameters: 



= 6/4l''^^i^2<3 / r^dr[-A,(r)A,(r)a,3(r) 



- {r)o^i^ {r)PH (r) " {r)Pi, {r)Pi, (r) 
-25i^{r)5u_{r)Si^{r) + {A^ {r)'yi2{r)5i^{r) + (5permutation)}], (18) 

where 

2 

71 



■yiir) 



k'dkP'J\k)gTi{k)ji{kr) 



&i{r) = \ I k'dkPl'\k)gTi{,k)ji{kr). 



(19) 
(20) 



The other type of NG which is sensitive to the bispectrum 
with a fold ed triangle configuration (l\ ~ I 2 — h /2) is also con- 
sidered bv lSenatore. Tassev. Zaldarriagal (|2009l '). This is called 
orthogonal-type NG and characterized by /j 



(ort). 



B.\°S3 -6/ir/.i.,,3y ^'rf^h3A,(r)A,(r)a,3(r) 

^Wh {r)m2 {r)Pi, (r) - 3a,, (r)/?,, {r)Pi, (r) 
-»5i,(r)5i,{r)5i,{r) -f 3{A, (r)7,,(r)5,3(r) + (5perm.)}]. (21) 

As d iscussed in Sec. 4.2 of ISenatore. Tassev. Zaldarriagal 
(|2009l l. we include the integration of r further than the 
last scattering surface to calculate the equilateral and 
orthogonal-type bispectra. This reduces the sensitivity of these 
types of NG. 

We also take into account the effect of unmasked point 
sources (e.g., radio galaxies) which generates an additional NG 
in observed CMB maps. Assuming them to be Poisson distribu- 
tion, the bispectrum has a constant value for all configurations 
of wavevectors: 



(23) 



-l^f^i + kU^(-)-^. 



v^2\^) = ^H,{v) + 



(24) 



K - 2SSi 
24 



-i (ki + is5„) Hs{u) - i [ku + iftSii) H,{v). (25) 



The kurtosis parameters are defined as 



K 



2{/|V/pVV)c + (lV/|^). 



^ _ ((vV)f)c 

= ota\ ' (26) 

where {...)c denotes the ensemble average of the connected 
part. 

For the local-type NG, the reduced trispectrum is written 

as 

T!^1^{L) = Ii,i.,Lli3UL I Y^NL J ridrir2dr2FL{ri,r2) 

XQii {ri)l3t2 {ri)ai3 {r2)l3t^ (r2) 
+3NL J r^drPi,{r)p,,ir)[ai,{r)Pi,{r) + /3i,{r)a,,ir)]Y (28) 

where 



-Fz,(ri,r2) = — / k dkP^{k)jL{kri)jL{kr2). 



(29) 



For the single-field inflation model with local- type NG (eq.[T]), 
tnl is equal to 36/jS,l'^'^/25. 

The kurtosis parameters are given by a summation over 
all configurations of trispectra. The full calculation of trispec- 
tra is, however, computationally very expensive and hence we 
estimate kurtosis values using Mon t e-Car lo integration of the 
flat-sky approximation. IMatsubaral (|2010l ) finds that the full- 
sky spectrum with its multipole configuration of {k} is well 
approximated by the fiat-sky spectrum with the wavelength 
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Figure 1. Three Skewness values (eq.[8]) for local-type f eg. 1151 ). equilateral-type feq. 1181 ). and orthogonal-type (eg. 1211 ) NGs and point 
sources as a function of the Gaussian smoothing scale 9. We add WMAP beam functions for V+W co-added maps and a pixel window 
function corresponding to each 8 listed in Table [l] 




9[arcmin] 0[arcmin] 9[arcmin] 0[arcmin] 



configuration of /i + 1/2: 



11112LI1314L 



(30) 



Tlie proportional factor Ii-^^i2Lli^iiL has non-zero value when 
both li -\- 12 + L and + I4, + L have even number. Note that 
the difference of 1/2 for all arguments of lengths I is very im- 
portant for the accurate estimation. The kurtosis in flat-sky 
approximation is given as 

xr(/i, I2M, h; I12, l23)Wih)W{l2)W{l3)W{h), 
where 



"0 



(31) 



2g2' 
A'lii 



^ it2 + 'iilili - ii2) 

32q^ 



(32) 



and the Gaussian window function in the flat approximation is 
now given by W{1) = exp[-{(/ + 1/2)6}^ /2]. The trispectrum 
T is symmetric against the arbitrary exchange of arguments 
among li, I2, I3 and I4 and then it is given as a sum of the 
reduced trispectrum T: 



T{li, h, I3, '4; '12, hs) — P{h, h, I3, U, '12) 

+P{li, I3, 12, li, hi) + P{li, I4, h, h, I23), 

and 



(33) 



P{li,l2,l:i,l4, L) = T{li, h, I3, 14, L) + T{l2, h, 13,14, L) 

+T{h, h, h, h, L) + Tih, h,l4, h, L). (34) 

As shown in Fig. [21 the quadrangle configuration is uniquely 
determined by five parameter spaces: three side lengths Zi, ^2, ^3 
and their open angles ^12, ^23. The other side length I4, two 
diagonal lengths l\2 and ^23 (see Fig. [2)) and /13 = |li + la] 
are written in terms of the five parameters. Fortunately, the 
configuration dependence of the trispectrum is smooth for the 
local-type NG model and hence Monte Carlo integration is 
applicable to estimate kurtosis values in shorter time. The in- 
tegration of the side lengths li, I2 and ^3 is done from 2 to imax 
given by Table [T] and that of ^12 and ^23 is done from to 2n. 
In equation H30p . the fractional values of the length {U} are 
allowed in the fiat-sky approximation of T, while the full-sky 



Figure 2. Trispectrum of two-dimensional CMB maps is calculated 
at each configuration of tetrahedron. A tetrahedron is uniquely de- 
termined by three side lengths Ii,l2,l3 and their open angles 9\2 
and 623. The other side length I4 and the diagonal lines I12 and I23 
used in the equation 1 133 II are also shown. 




Table 1. Resolution scale Afgjjc for HEALPix (total number of pix- 
els is given by 12N^^^^) s 
at each smoothing scale 1 



els is given by UN^.^J and the maximuni values of multipole /n 



9 [arcmin] 




^max 


100 


128 


136 


70 


128 


196 


40 


256 


340 


20 


512 


684 


10 


512 


750 


7 


512 


1024 


5 


512 


1250 



trispectrum T'^i^ (L) is only given for a set of integer values k . 
We use the round-off values of {h} to get T. We get the rea- 
sonable values of kurtosis parameters by 0(10^) calculation, 
which takes much shorter time than the full calculation of the 
trispectrum. 

For the trispectrum, we also consider the effect of point 
source as 



where t^'^^'' is a constant value. 



(35) 
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3 COMPARISON OF PERTURBATION THEORY 
WITH NUMERICAL SIMULATIONS 

3.1 Computation of Minkowski Functionals from 
CMB maps 

Computation method of the MPs of CMB maps are d e scribe d 
in Appendix A.l. of iHikage. Komatsu fc Matsubaral (|2006l ). 
The range of u is from —3.6 to 3.6 and the binning number 
is set to be 36 for each MP with the equal binning width 
Ai^ = 0.2. The numerical estimations of 2nd and 3rd MPs 
do not completely agree with th e anal ytical predictions even 
in Gaussian fields. iLim fc Simoiil (|2012l ) found that the devia- 
tion comes from the approximation of the delta function with 
a finite difference. In Gaussian fields, the correction terms of 
the finite binning effect of the 2nd and 3rd MPs are given by 



1 



dxVk{x) 



(36) 



where Av is the binning width of u. We subtract the correction 
terms from the measured MPs. 



3.2 Comparison of Perturbation Theory with 
Non-Gaussian CMB Maps 

Actual CMB maps have various observational effects such as 
survey geometry, inhomogeneous noise, which may cause sys- 
tematic uncertainty in the NG measurements from MPs. We 
take into account such observational systematics by construct- 
ing realistic CMB simulation maps and test if the perturbation 
works for actual observations. 



3.2.1 Full N on- Gaussian Simulations 

We employ 1000 realizations of si mulated CM B maps with 
a local- type NG available in public (lElsner fc Wa ndclt 2009f ). 
These simulation maps include full radiat ion transfer function . 
We make mock CMB maps as described in lHikage et al.l (|2008l ): 
making Gaussian CMB maps with their input power spectra 
following WMAP 7-year cosmology and including the beam 
function for each differencing assembly (DA), we add Gaussian- 
distributed noise to each pixel with the standard deviation 
of a/Nohs in WMAP 7- year observations. We co-add noise- 
included DA maps with the inverse weight of averaged noise 
variance and then mask the area outsi de the KQ75 ma sk of 
Galactic foregrounds and point sources (jCold et al.ll201l1 ). 

Pig. [3] (left-hand panels) shows a comparison of two vari- 
ances (eq.U), three skewness parameters (eq.[S]), and four kur- 
tosis parameters (eq.[5S]) with f^'^'^^ = 300 between analytical 
estimations and simulation results at different values of 6s. 
The simulations do not include the pnl term and thereby the 
second-order NG comes from the square of the /nl term, or 
tnl. The error-bars represent the la dispersion of simulation 
results divided by the square root of 1000, that is, the number 
of realizations. The right-hand panels show that the 2nd-order 
correction of each MP at S = 10 arcmin. In the plot of MPs, 
we subtract Gaussian and Ist-ord er perturbative corr ection to 
focus on the 2nd-order correction (jlfikage et al.l (|2008l ) have al- 
ready shown that the Ist-order perturbative correction of MPs 



find that the 2nd-order perturbative corrections of MPs also 
agree with simulation results very well even including the ob- 
servational effects. 



3.2.2 Sachs-Wolfe approximations 

We also compare the 2nd-order correction due to (;nl type us- 
ing NG simulation maps in Sachs- Wolfe approximation where 
the NG is locally given by 



ATb 
T 



3/, 



(loc) 



( ATg 
\ T J 



V 



(37) 



The angular power spectrum is set to be l{l -f 1)C; /27r = 
at / <= Zmax where Zmax is given in Table [T] at each 0. We 



IQ- 



generate 6000 realizations with f^^'^'' ~ 100 and (;nl ~ 10^ . Por 
these values, the Ist-order correction from f^^'^^ and the 2nd- 
order correction from pnl have comparable amplitudes, but the 
contribution of tnl ~ /nl'^'^ is negligible. Pig.|3]shows a similar 
plot to that of Pig. [3] but for the comparison with simulations 
in the Sachs- Wolfe approximation. The right-hand panel shows 
the MPs subtrac ting only Gaussian terms at 6* = 10 arcmin. 
iMatsubaral (|2010l ') showed that the perturbation formulae work 
very well in the Sachs- Wolfe approximation. Here we find that 
the perturbation also works even including the various obser- 
vational effects such as survey mask, inhomogeneous noise and 
beam window function. The excellent agreement indicates that 
the flat-sky approximation and Monte Carlo integration also 
work very well at a wide range of scales. 



3.3 Distribution of NG parameters estimated from 
simulations 

We estimate the bestfit values and the errors of each NG pa- 
rameter by the least chi-square fitting of observed MPs with 
the perturbative predictions: 



X 



E 



^^(obs) _ ^^(thcory) ^ j 1 j^(obB) _ ^(theory) ^ ^gg^ 



due to f^l'^^ agree with the results from NG simulations). We 



where pk denote NG parameters such as /nls, tnl and (?nl 
and subscript i and j of MPs denote different bins of threshold 
V, different kinds of MPs, and different smoothing scales 6s. 
We estimate the covariance matrix of MPs using 6000 realiza- 
tions of Gaussian CMB maps including WMAP7 observational 
effects described in the previous section. Por the chi-square 
measurement, we reduce the binning number to 18 from 36 
original bins. The reduction of binning does not affect the re- 
sults, which means that the binning number of 18 is enough 
that the result is converged. The maximum number of the side 
length of the covariance matrix is 270 corresponding to 18 bins 
of X three kinds of MPs x five different 6. We have checked 
that 6000 realizations are enough for the fitting results to be 
converged. 

Assuming that the covariance derivative ddj/dpi is neg- 
ligible, the Pisher matrix is simply given by 

^ dpk dpk' 

We use NG CMB simulation maps and test if the above chi- 
square estimations generate the expected distribution of k- 
th NG parameter with the mean of the input value and the 
error expected from the Pisher matrix as [{F-')kk]''^. Pig. 
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Figure 3. Left: varianc e, skewncss and kurtosis of NG simulations with /^'l'^' = 300 using CMB simulation maps with local-type NG 
l lElsner fc Wandeltll2009l) (symbols) . The lines show the theoretical predictions based on the perturbation theory. Right: second-order correction 
of MPs in the same NG simulation maps smoothed at 5 = 10 arcmin (symbols) by subtracting the Gaussian term V, and the first-order 



correction V, 



(1) 



, ^ . The NG simulations do not include the component and thereby the second-order correction comes from t^l corresponding 

1.296 X 10^. For comparison, perturbative predictions are written in the lines. The simulations include the WMAP beam 
and noise for the V-|-W co-added map and a pixel window function at each 9. 
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[5] shows the distribution of best-fit values of /nl , tnl from 
full NG simulations and pnl from NG simulations in Sachs- 
Wolfe approximation. The details of NG simulations are writ- 
ten in the previous subsection. Averaged values are respectively 
/nl = 100±47(-100±48) against the input of 100(-100)±48, 
and tnl/IO" = 12.3 ± 8.6(12.4 ± 9.1) against the input of 
13±8.5, and 5nl/10^ = 10±1.76 against the input of 10±1.5. 
We find that our method well reproduces the input values of 
NG parameters and the error estimations using the Fisher 
matrix (eq.[39]) provide reasonable measurement error. The 
mean value of tnl is found to be underestimated by 5%, which 
may be due to the incompleteness of the theoretical estimation 
of the kurtosis parameters. The systematic error is, however, 
much smaller than the statistical error of tnl and the effect on 
the final result is small. 



parameters carry the NG information and hence their mea- 
sured values can be directly used to limit the NG parameters. 
For consistency checks, we compare the constraints on NG pa- 
rameters estimated from MFs with those estimated directly 
from skewness and kurtosis parameters (hereafter we call them 
"moments") using WMAP mock simulation maps. Fig.|6]shows 
the distribution of f^'^^ estimated from MFs (red), and those 
from skewness values (blue). We find that the sample variance 
are comparable and thus MFs and the moments have similar 

We also plot the difference of the 



power to constrain f^"^^ 



bestfit values of f'^'^'^ from the MFs and the skewness divided 
\/2. The limits from the MFs and the moments do not com- 
pletely agree because the weights on skewness parameters in 
MFs are not equal. The difference is smaller than the dispersion 
of the distribution of /nl , which means that the measurements 
of MFs and moments are strongly correlated. 



3.4 Comparison of constraints from MFs and 
skewness /kurtosis 

Perturbative corrections of MFs are determined by skewness 
(eq:[8]) and kurtosis parameters (eq.[26]) upto 2nd order. The 
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Figure 4. Left: same as Fig.[3]but for NG simulations in the Sachs- Wolfe approximation. The NG parameters are /^'l'^' = 100 and ^nl = 10^ 
in which the contribution of t^l is negligible. Right: comparison of the NG correction of each MF between the perturbative theory (solid 
lines) and the NG simulations (filled circles). First- and second-order perturbative corrections are shown with the dotted and dashed lines, 
respectively. 
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Figure 5. Left: distribution of the bestfit values of f^'^'^'' estimated from 1000 NG simulation maps in which the input values of f^'^'^^ = ±100; 
Middle: distribution of t^l estimated from 1000 NG simulation maps with the input values of /^'l'^' = ±300, that is, the corresponding 
tnl = 1.296 X 10^; Right: distribution of jnl estimated from 1000 NG simulations in the Sachs- Wolfe approximations. The input values of 
NG parameters are /^'l'^' = 100 and gNL = 10® . The lines show the Gaussian distribution with its mean of the input values and the dispersion 
given by the square root of the inverse of the Fisher matrix ■ 
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Figure 6. Distribution of the bestfit values of /j 
input value /JJ^^ = 100 using 1000 NG simulation maps. We obtain 
the bestfit values using two different measurements: MFs (red lines) 
and three skewness parameters (blue lines). The differences between 
these two measurements are also plotted with the yellow lines. 
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4 APPLICATION TO WMAP 7- YEAR DATA 

In the previous section, we show that the perturbation works 
even including the various observational effects. We apply the 
perturbative formulae to WMAP 7-year temperature maps and 
give limits on the NG parameters. Fig. [7] shows the MFs of 
WMAP 7-year V+W co-added maps at different smoothing 
scales. We subtract Gaussian term Vj, from the observed 
MFs to focus on the deviation from Gaussian. For comparison, 
we plot the perturbative formulae with the bestfit values of 



/nl'^ '^NL, and ^nl- The Ist-order {f^'^ ) and 2nd-order (tnl 
and (/nl) contributions are shown with dotted and dashed Hne, 
respectively. The differences are consistent with zero for all 
MFs at different smoothing scales. Fig. [8] shows the skewness 
and kurtosis parameters of the same WMAP data at different 
smoothing scales, which are also consistent with zero. 

Table [2] lists the limits on /nls, tnl and gNL. The con- 
straining power is strongest at 6'=10 or 7 arcmin scales be- 
cause the noise is dominated at smaller scales. Maps of differ- 
ent smoothing scales have different scale information and hence 
combining results from different 6 maps provide stronger con- 
straint. These constraints take into account the point source 
effect and its contribution is marginalized over. All types of pri- 
mordial NG parameters we consider are consistent with zero. 
The constraints from MFs ar e weaker than the opti mal estima- 
tor based on the bispectrum |Komatsu et al.|[201ll ) and on the 
trispectrum (|Smidt et aI.ll201Cll ). This is because that the skew- 
ness and kurtosis parameters lose configuration information on 
the bispectrum and the trispectrum. Again we stress that the 
consistency check from MFs are important to check the sys- 
tematic effects. It may be interesting that the bestfit value of 
Tnl from MFs is slightly inclined to be negative because the 
negative value of tnl is not allowed for all of multi-field in- 
fiation models by the inequality relation tnl > 36/25/4l'^'^ 



Figure 8. Three skewness and four kurtosis parameters measured 
from WMAP 7-year data at different smoothing scales of 6. The 
definitions of the skewness and kurtosis parameters are given in the 
equations ^ and I I26I I. respectively. 
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ama fc Yamaguchill2008l ). This inclination is stronger be- 
fore marginalization of point source effect as seen in Table |3] 
however the significancy is still very small. We also give the 
limits from different frequency bands listed in Table|4] The dif- 
ferences between different frequency maps are within 1 sigma 
of statistical error, which means that the frequency-dependent 
systematics such as Galactic foreground do not affect our re- 
sults so much. The results from MFs and moments (skewness 
and kurtosis) are also consistent and their difference is within 
1(T statistical error. 



5 SUMMARY AND DISCUSSIONS 

We first apply the perturbative formulae of MFs including 
second-order NG to WMAP 7-year data and give limits on tnl 
and (jNL as well as local-type, equilateral-type and orthogonal- 
type /nl. Consistency check using different estimators are im- 
portant to obtain more robust results because different es- 
timators are sensitive to different aspects and systematics. 
We find no evidence of NG from any type of NG compo- 
nents and then obtain the limits on each NG parameter: 
/^'r' = 20 ± 42, /(jI^' = -121 ± 208, and f^Z'" = -129 ± 171, 
tnl/IO"* = -7.6 ± 8.7 and c/nl/IO^ = -1.9 ± 6.4. Our result 
is consistent with the previous works using the estimators of 
bispectrum and trispectrum. Constraining tnl is important 
for testing multiple infiation models which must satisfy the 



inequality condition tnl > 36/25/nl. Our limit on tnl is con- 
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Figure 7. Three MFs for WMAP 7-year temperature maps at different 6=40, 20, 10, 7, and 5 arcmin from the top to bottom. The Gaussian 
term V^. is subtracted to focus on the deviation from Gaussian. For reference, the first and second-order perturbative corrections with the 
bestfit values of /^'l'^' , tnl a^nd 9nl a-re plotted, respectively. 
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sistent with zero, but the bestfit value of tnl is inclined to be 
negative even after subtracting point source effects. Upcoming 
CMB experiments such as Planck will give statistically better 
results. 

Large-scale structure offers another test to study the 
presence of primordial NG. Scale-dependent bias of halo 
clustering al s o pro vides a powerful probe of primordial NG 
jPalal et al.l |2008| ) and observational constraints —29 < 
fi'r' < 70 (95%C L.) are obtained jSlosar et all l2008h . 
iDesiacgues fc SeliakI (I2OI0I ) give constraints on (;nl as —3.5 x 
lO'^ < c/NL < 8.2 X lO'^ (95% C.L.) from the halo mass function 
and halo bias. Combining galaxy bispectra is also useful to 
constrain the primordial NG including the higher-order NG s 
Tnl and (7nl (|Jeong fc Komats"ull2009l : [Nishimichi et al.ir2010l ). 

The secondary effects like a coupling between the in- 
tegrated Sachs-Wolf e effe ct and the gravitational lensing 
(|Goldberg fc Spergell 1 19991 ) may contaminate the measure- 
ment; however, the est imation is small of the order of /nl ~ 3 
(|Komatsu et al.l I2OIII ). In this paper we give limits on the 
three types of /nls and tnl and (7nl locally given. In gen- 
eral, there is still a wide range of NG such as isocurvature 



NG (|Kawasaki et aklbOOSlTHikage et al.ll2009l). The equilatera l 
type of (7NL has been also given by Mizuno fc Kovainal (|2010l ). 
The application to the other types of NG is a future work. 
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Table 4. Limits on the NG parameters in different frequency bands. All of the limits are obtained from MPs (Upper) or Moments (Lower) 
by combining different ds maps and the point source effect is marginalized over. 
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